Vanishing Spin-Hall Conductivity in 2D disordered Rashba electron gas 
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For the two-dimensional ideal electron gas with the Rashba spin-orbit interaction and in the 
presence of non-magnetic short-ranged potential impurities the spin-Hall conductivity aaH is found 
by direct microscopic calculation. Within the semiclassical approximation h/r <g ep and A <C ef 
the value of Gsh ts zero for arbitrary ratio of the spin-orbit splitting A and the inverse elastic 
scattering time 1/r. 



PACS numbers: 72.25.-b, 72.10.-d, 72.15.Gd, 73.50.Jt 

Recently it has been proposed^ that a dissipationless 
spin current can be generated in response to an elec- 
tric field in semiconductors with the spin-orbital inter- 
action. For the case of an ideal two-dimensional (2D) 
electron gas with the Rashba coupling, Sinova et ali^ 
have found a spin-Hall current of the transverse (2) spin 
component as a response to an in-plane electric field E^, 

— UsH^nuE^, with the "universal" (independent upon 
the spin-orbital band splitting A) spin-Hall conductivity 
CTsif = e/87r?l. In the presence of disorder and for the 
simplest standard model it was demonstrated a cancella- 
tion of the spin-Hall effect even in the case of arbitrary 
week disorder In this communication we provide a 

fully microscopic calculation of the spin-Hall conductiv- 
ity for the generalized model of non-parabolic spectrum 
and arbitrary momentum dependence of the Rashba ve- 
locity iQ}, within the semiclassical theory of disordered 
conductors, i.e. in the limit U/t <C ep^ A <C ep- We 
show, that UsH = independently upon the parameter 
tA/Tl sim ilarly to the case of the parabolic spectrum 
Ref. MM 

2D isotropic Rashba gas is an electron system with the 
broken inversion symmetry. In this case an electric field 
perpendicular to the plane could arise. It has no effect on 
the electron orbital motion but it couples to the electron 
spin via a relativistic spin-orbit interaction known as the 
Rashba term. The Hamiltonian of an electron consists of 
the kinetic energy term and the Rashba term: 



(1) 



where ol{p) is the Rashba velocity, u' {i — x, y, z) are the 
Pauli matrices and a, (3 are the spin indices. The Hamil- 
tonian can be diagonalized by the unitary matrix: 



U{p) 



V2 \ 



(2) 



where ifp is the angle between the momentum p of the 
electron and the x-axis, with the eigenvalues: 



(3) 



The eigenvalues A = ±1 of the chirality operator and 
p constitute the quantum numbers of an electron state 



{p, A) . Fermi circles of the Rashba gas with the different 
chiralities are split: pf± = Pf{^^ Oi{pF)/v{PF)), where 
Fermi momentum pp solves the equation: £{pf) = /x, 
where /i is the chemical potential; v{p) = de{p)/dp is the 
band velocity of the electron. The electron velocity in 
the chiral state is de\{p)/dp. We assume a{pp) <^ v{pp), 
and neglect corrections of the order a/ v. The spin-orbital 
splitting is then: A = 2ppa{pp). The density of states 
on the two Fermi circles differs as: i>± — i^_f(1 ± ap/vp), 
where vp — pp /2'Kv{pp). Contrary to the case of the 
parabolic spectrum and Rashba velocity independent on 
the momentum, for the generalized model the Fermi 
velocities are different on the two Fermi circles: vp^ — 



Vp 



2ap(j^^ 1) 

use the units where h = 1. 

We consider the 2D ideal (non-interacting electrons) 
Rashba electron gas with the Hamiltonian: 



2pp 



In the following we 



Hr ^^al{p) hap{p) ap{p), 



(4) 



at zero temperature, al^ip) and ap{p) are the electron 
creation and annihilation operators. Electromagnetic 
vector potential A couples to the orbital motion of the 
electron according to the transformation: p^ p — eA/c, 
in the Hamiltonian ^ . Variation of the Hamiltonian Q) 
with respect to A gives the electric current operator: 
— J2pO'l,ip){ju)ai3ai3{p), where the one-particle cur- 
rent operator reads (it is actually a velocity — ei)^): 



. ^P<.^ , d{pf,a{p)) 

v[p)—dafi H 5 

P dpy 



(5) 



with V = x,y being the spatial index and e^'^ is the 3D 
totally antisymmetric tensor. 

Under the non-uniform SU(2) electron spinor transfor- 
mation: aa{r) ^ Uai3(r)ai3{r), the Hamiltonian |0J be- 
comes dependent on the SU(2) "spin electromagnetic" 
vector potential A^^ — A^fi^ -I- A^^cr^ , where A° coin- 
cides with the physical electromagnetic potential and 
= -iTr(cr*f/+a^?7)/2. Although this latter poten- 
tial is a pure gauge and has no physical consequences, 
variation of the Hamiltonian |0J with respect to it de- 
fines the spin current of i-component of the spin along 



2 



the direction /i: = Y.p'A^{'P){fti)al30'p{p), where the 
one-particle spin current operator reads: 



(6) 



Our definition of the spin current lO coincides with the 
definition followed in Ref. 0: = {'I'fi'^^ + o'*'^^)/2, but 
differs from the definition followed in RefsiSi^ by a factor 
2, which makes our value of the spin-Hall conductivity 
being twice as big as that in the literature liiSi^ 

The interaction of electron with short-ranged non- 
magnetic impurities at positions numerated by the 
index i, is described by the impurity Hamiltonian: 



(7) 



where u{r) is a short-range impurity potential. We as- 
sume it to be sufficiently weak in order for the Born ap- 
proximation to be valid. In this limit the impurity model 
([7|) is equivalent to the model of the Gaussian random 
potential. We expand the electron Green function aver- 
aged over the realizations of the disorder potential pur- 
turbatively in power of the Hamiltonian (jjj) using the 
diagrammatic procedural. It is a sum of diagrams where 
a chain of electron bare Green functions is separated by 
impurity "crosses" . Two "crosses" are connected by the 



averaged impurity line: riy 



l/27ri^T, where 



the density of impurities, and r is the scattering mean 
free time. A "cross" does not change the electron spin 
and the electron frequency since the electron scattering 
off impurity is elastic. Therefore the impurity line carries 
zero frequency. Diagrams with crossings of two or more 
impurity lines are small as powers of the ratio 1 /epT 1. 
The averaged Green function is a two by two matrix in 
the spin space and it is a solution of the Dyson equation: 



V 



pi 

(8) 

It can be conveniently transformed into the chiral ba- 
sis by the unitary matrix lJ(j>) |(2J). The retarded and 
the advanced averaged Green functions are diagonal in 
the chiral basis: G^^^\t,p) = G^^'^^ 
solution to the Dyson Eq. (jSJ reads^: 



G\'^''^'{e,p)5x'\, and the 



1 



ex{p) + ^Ji + i/2T 



(9) 



where r is explicitely independent of the chirality. The 
advanced Green function is a complex conjugate of the re- 
tarded one: G^{e,p) = {Gf{e,p)}*. The Green function 
in the spin basis is a non-diagonal two by two matrix: 



G\f''\e,p) G\f^\e,p) 
G\f^\e,p) G[f^\e,p) 



(10) 



where (omitting for a moment the frequency and momen- 
tum notations): 



TT 

n 

{R,A) 
U 



G 



(Gf-^) 



G 



(R.A) 



)/2, 



(11) 



with the chiral Gj^''^'' being defined in Eq. 

In order to calculate the current, induced in the elec- 
tron system by an electric field, we use the Keldysh 
techniqueiS. Our result given by the Eq. (|14|l is well 
known but we derive it here for consistency. The aver- 
aged Keldysh Green function is a four by four matrix 
Q (p, e) that can be conveniently factorized into a two by 
two Keldysh matrix whose elements are matrices in the 
spin space themselves: 





G- 












-{ 



l-N{() -N{e) 
l-N{e) -N{e) 



G^(p,e) 



+ 



N{e) N{e) 
-1 + N{e) -1 + N{e) 



G^(p,e), 



(12) 



where the electron distribution N{e) is proportional to 
the unit matrix in the spin space. 

We choose the gauge for the uniform electric field 
E{t) — E{VL)e~^^^ to be a time dependent vector poten- 
tial A{t) = A{Q)e-'^\ where A{Vt) = -icE{n)/n. Us- 
ing the Keldysh techniqueifi we average the spin current 
operator over the electron state perturbed by the elec- 
tromagnetic Hamiltonian: Hem = ~- J cP^ju{T^)^u{t), 
in the first order of the perturbation theory. The spin- 
Hall conductivity asH is then found from the relation- 
ship: (j^(f2)) = e^^CT^//(f7)£'i.(fJ), as: 



—y 



27r 



-Tr 



-+ 



(13) 

where is the four by four matrix given by the direct 
product of the Pauli matrix in the Keldysh space and 
the unit matrix in the spin space, the current operators 
in Eq. (|13|) are the direct product of matrices ^ and 
the unitary matrix in the Keldysh space. Tr in Eq. H13() 
operates only in the spin space whereas the indices of — h 
element corresponds to the Keldysh space. Substituting 
the Green function Eq. (|12|) . we obtain: 

p 

X (G«(e + r!, p) - G-4 (e + r!, p)) {p)G^ (e, p) 
+ G^(e + n,p)Up)N{e) {G%,P) - G^{e,p)) )] ),(14) 

where the brackets indicate averaging over the disorder. 
In the non-crossing approximation the average in Eq. H14|) 
is given by the sum of the one-loop and the ladder dia- 
grams shown in Figs.m^l 



FIG. 1: The spin-Hall conductivity is given by a one- loop 
diagram. 



FIG. 2; The vertex correction to the spin-Hall conductivity 
is given by a sum of non-crossing ladder diagrams. 



First, we calculate the one-loop diagram Fig. ^ and 
denote this part of the spin-Hall conductivity as cr^//. 
It corresponds to the Eq. H14(l with all Green functions 
being substituted by the averaged Green functions 
The second line in Eq. (|14|l contains the imaginary part 
of the Green function G^{€ + Q,p) and the corresponding 
integral is convergent, therefore we change e to e — f2 in 
this second line. At T = the Fermi-Dirac distribution 
function reads: N{e) — d{—£)- We take the integral over 
e and in the zero-frequency limit O ^ we find: 



1 Hp) 

p 

1 



1 



2pa{p) 



is{C) - s{0) 

f 1/4t2 



2r (C2 + 1/4t2)(^2_|_ 1/4^2) 



(15) 



where S{x) = arctan(a;/2T), ~ e{p) — pa{p) — /i and 
S. — e{p) +pa{p) — /i. In the large volume limit we substi- 
tute y ^ - — > / ' ^^'^ then evaluate the integral over 
p in H15|l in the limit of large Fermi circle 3> A, l/r. 
This procedure, known as the semiclassical approxima- 
tion, expresses the momentum p in terms of the quasi- 
particle energy ^ : 



d^p 

(27r)2 " 271 v{p) J_ 



il + R{p)Od^ 



dip 
2^' 



(16) 



where Rip) ~ p ^(p) ~ m ^{p)v ^(p). The result 
reads: 



'sH 



in ^ 



1 



1 + (At) 



(17) 



with the first and the second terms corresponding to the 
first and the second terms of Eq. lfT31) , respectively. 

An important observation is that the result (|17|l coin- 
cides exactly with the result obtained from those terms 
in Eq. (|14|l that contains one retarded and one advanced 
Green functions. These two terms are proportional to 
dN{e)/de ~ — <5(e) and all integrals are explicitely con- 
fined to the vicinity of the Fermi circle. Therefore the 
spin-Hall conductivity unlike the usual Hall conductivity 
is determined by the quasi-particles around the Fermi 
circle and not by the entire Fermi disk. 

The ladder diagrams shown in Fig. |21 represent the ver- 
tex corrections to the current. Additional impurity lines 
improve the convergence of the integral in Eq. (|14|l . As a 
consequence, vertex corrections to the terms in Eq. H14|l 
with the two advanced or with the two retarded Green 



functions vanish as max(l/T, A)/eF ^ 1. Therefore we 
consider only the vertex corrections to the terms with 
one advanced and one retarded Green's functions as it 
is shown in Fig. 13 For these diagrams the semiclassical 
approximation ]l<i\i is valid. 

The sum of ladder diagrams with n = 1, 2, ... impurity 
lines is given by the expression: 



„lad 

f^sH = - 



d^p 

{2nf 



Tr 



4G«(0,p)j,(p)G^(0,p) , (18) 



where the sum of n = I..00 vertex corrections to the 
current jy (p) (with at least one impurity line) is denoted 
by the matrix Jy. In the spin basis and for short-ranged 
impurity potentials it does not depend on the electron 
momentum p and satisfies the transfer matrix equation: 



Jy 



1 



d^p 

2tttv J (27r)2 



(19) 

where the Green functions G'^^'^^ are given by Eq. ll()|l . 
The "full" spin current operator with all vertex cor- 
rections included: j^ip) + Jy , is represented diagram- 
matically in Fig. |3l In the equations for the cur- 
rent operators (jS] IHJl we expand the electron veloc- 
ity v(j>) = v{pp) + ^/ {v{pF)'m{pp)), where m~^{p) = 
dv{p)/dp, to the first order in the deviation from the 
Fermi circle: ^//x, small in the semiclassical approxima- 
tion. We also expand the spin-orbital splitting: pa{p) = 
a{pp) (pp + [1 + (j)p /aF)da/dpF)^/vF), in the Green 
functions. Wc then evaluate Eq. (|19|l in the semiclassical 
approximation H16|l neglecting odd powers of ^: 



^ {^^"^' (^.^),, + [2 + (ATn (^.^)^}^' 

[Jl)^^ - |-z.(pF)Ar+ [2 + (Ar)^] 

[Jl)^^ = L{pF)AT+[2 + (Arr] {J^)\b, (20) 



where B ~ 5 [l + (At)^] ^ . From the first two lines 
of Eqs. we find: (J^)tt — {Jy)ui whereas from 

the last two lines wc find: (J^)xt = ^v{pf)/{At) and 
(J^)lj^ = -iv{pp)/{AT). The integrand of Eq. lO does 
not depend on the diagonal elements of the matrix Jy and 
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FIG. 3: The vertex of the spin current with the vertex cor- 
rections taken into account: jyip) + Jy. 

therefore we set them to zero: Jy = cr^ v{pf)/ ^t. As it 
was expected the vertex corrections are proportional to 
the scattering rate. Integrating Eq. (|18|l in the semiclas- 
sical approximation H16|l over ^, we finally find the ladder 
part of the spin Hall conductivity (Fig. : 

""^i(-' + TT(W)^ 

Remarkably all derivatives of a{p) and v{p) over p have 
canceled out. 

The spin-Hall conductivity is the sum of Eqs. ifTTll^ 
and is zero: 

asH = a% + a't^ = 0. (22) 

It explicitely does not depend on the impurity scattering 
time T. But, we observe a discontinuity between the spin- 
Hall conductivity in the clean system asH = e/Sirh and 
the spin-Hall conductivity Eq. H22() in the presence of the 
infinitely small amount of non-magnetic scatterers. As it 
was shown in Ref. 6, this discontinuity is related to the 
dissipation in the system, which gives rise to the dissi- 
pative part in the spin-Hall conductivity crf^ = ^e/8nh, 
which cancels the reactive part = e/SirTi. 

In an analogous calculation, for the generalized 
model (0, we find the magneto-electric effect - magneti- 
zation induced by the electric field: 

{Sn = e^'^T^i?., (23) 
Zttvf 

where S'^ = 1/2 J2p^ai'P) ^afs^f^ip) the total spin op- 
erator and the total magnetization {S'^) is in agreement 



with Ref. H IIII12I One should notice that the steady 
in-plane magnetization (|23|l is a consequence of the zero 
spin-Hall effect (1^ . 

Non-zero spin Hall conductivity would result in a non- 
steady in-plane magnetization. This follows from the evo- 
lution equation of the total spin of the system Si, and the 
commutation relationship: 

-^p^^[H,S'^]^^aip)-^J^^, (24) 



where fi = x,y; is the total spin current operator, de- 
fined in lO . Remarkably, Eq. H24I) is valid for an electron 
system with any non-magnetic disorder and any electron- 
electron interaction. Moreover, valid is a more general 
equation for the local spin density evolution: 

^is'^ir) + l^J^ + am [e^y^g - e^^^^y) = 0, (25) 

where ^ — x,y, z; in derivation of Eq. (|25() we considered 
the standard model of parabolic spectrum. If the last 
term were zero then the above Eq. (|25|l would be the spin 
conservation equation. But in the Rashba metal neither 
component of the spin is conserved. 

To conclude, we have extended the result of Inoue et 
al. Ref.i cTsH = for the case of the arbitrary electron 
dispersion, arbitrary strength of disorder and arbitrary 
momentum dependence of the Rashba velocity a{p). Our 
result agrees with Ref. 
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